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Oh 1 ABSTRACT 

i 

2 . The VISTA near infrared survey of the Magellanic System (VMC) will provide deep YJK S photometry reaching stars in the oldest turn-off 
point all over the Magellanic Clouds (MCs). As part of the preparation for the survey, we aim to access the accuracy in the Star Formation 
^ History (SFH) that can be expected from VMC data, in particular for the Large Magellanic Cloud (LMC). To this aim, we first simulate 
1 1 , VMC images containing not only the LMC stellar populations but also the foreground Milky Way (MW) stars and background galaxies. The 
simulations cover the whole range of density of LMC field stars. We then perform aperture photometry over these simulated images, access the 
expected levels of photometric errors and incompleteness, and apply the classical technique of SFH-recovery based on the reconstruction of 
colour-magnitude diagrams (CMD) via the minimization of a chi-squared-like statistics. We verify that the foreground MW stars are accurately 
recovered by the minimization algorithms, whereas the background galaxies can be largely eliminated from the CMD analysis due to their 
particular colours and morphologies. We then evaluate the expected errors in the recovered star formation rate as a function of stellar age, 
■ SFR(r), starting from models with a known Age-Metallicity Relation (AMR). It turns out that, for a given sky area, the random errors for ages 
older than ~ 0.4 Gyr seem to be independent of the crowding; this can be explained by a counterbalancing effect between the loss of stars 
due to a decrease in the completeness, and the gain of stars due to an increase in the stellar density. For a spatial resolution of ~ 0.1 deg 2 , the 
random errors in SFR(f) will be below 20% for this wide range of ages. On the other hand, due to the smaller stellar statistics for stars younger 
^} [ than ~ 0.4 Gyr, the outer LMC regions will require larger areas to achieve the same level of accuracy in the SFR(f). If we consider the AMR as 
unknown, the SFH-recovery algorithm is able to accurately recover the input AMR, at the price of an increase of random errors in the SFR(f) 
by a factor of about 2.5. Experiments of SFH-recovery performed for varying distance modulus and reddening indicate that these parameters 
K> can be determined with (relative) accuracies of A(ra-M)o ~ 0.02 mag and AEb-v ~ 0.01 mag, for each individual field over the LMC. The 
j_J ■ propagation of these latter errors in the SFR(/) implies systematic errors below 30%. This level of accuracy in the SFR(/) can reveal important 
' imprints in the dynamical evolution of this unique and nearby stellar system, as well as possible signatures of the past interaction between the 
' MCs and the MW. 

Key words. Magellanic Clouds - Galaxies: evolution - Surveys - Infrared: stars - Hertzsprung-Russel (HR) and C-M diagrams - Methods: 
numerical 
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1. Introduction 

Determining the star formation histories (SFH) of the 
Magellanic Clouds (MC) is one of the most obvious goals in the 
study of nearby galaxies, for a series of reasons. First, this SFH 
does likely keep record of the past interactions between both 
Clouds and the Milky Way d01senlll999l:lHoltzman et alJll999t 
ISmecker-Hane et ail I2002I : Harris & Zaritskvl I2004 ), whic h 
are still to be properly unveiled dKallivavalil et al. l2006bllaL 



Besla et al.l2007l:IPiatek et alj2008l) . Detailed SFH studies may 
also provide unvaluable hints on how star formation is trig- 
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gered and proceeds in time, from the smallest to galactic-size 
scal es, and how these processes depend on dynamical effects 
(e.g. lHarris & Zaritsky||2"007l:lHarrisll2007ah . 

The Magellanic Clouds are also a rich laboratory for stud- 
ies of star formation and evolution, and the calibration of pri- 
mary standard candles, thanks to the simultaneous presence 
of a wide variety of interesting objects such as red clump gi- 
ants, Cepheids, RR Lyrae, long period variables, carbon stars, 
planetary nebulae, the tip of the red giant branch (RGB), dust- 
enshrouded giants, pre-main sequence stars, etc. Although the 
system contains several hundreds of star clusters for which age 
and metallicity can be measured, the bulk of the interesting stel- 
lar objects are actually in the field and irremediably mixed by 
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the complex SFH, and also partially hidden by the presence 
of variable and patchy extinction across the MCs. Unveiling 
this complex SFH may help in calibrating stellar properties - 
like luminosities, lifetimes, periods, chemical types, etc. - as a 
function of age and metallicity. 

In the last two decades, many authors demonstrated 
that recovering the SFH of the MC from optical photom- 
etry is indeed feasible and well worth of the effort. Such 
works are, usually, based either on deep Hubble Space 
Telescope (HST) photometry reaching the oldest main se 
quenc e turn-off for small MC areas 



[e.g^ Gallagh er et al,. 

1996 1 iHoltzman et all fl999l lOlsenl 1 19991: lElson et all Il997t 
Smecker-Hane et alj|2002t lArdeberg et alJll997HDolphin et al] 
Javiel et al. Il2005h . or on relatively shallo w ground-based 



2001 



photometry covering larger areas over the M Cs ( Stappers et al . 1 



1997 



2001 



Gardi ner & Hatzidimitrioul 1 1992c lHa rris & Zaritsky 
2004 . Only in very few cases (e.g. iGallart et all 12004 



Noel et al 2007)' have the ground-based optical photometry 



been deep enough to reach the oldest main sequence turn-offs. 

The VISTA S urvey of the Magellanic System^ (VMC, see 
Cioni et al.ll2007[ Cioni et al., in preparation) is an ESO pub- 
lic survey project which will provide, in the next 5 years, crit- 
ical near-infrared data aimed - among other goals - to im- 
prove upon present-day SFH determinations. This will hope- 
fully open the way to a more complete understanding of how 
star formation relates to the dynamical processes under way in 
the system, and to a more accurate calibration of stellar mod- 
els and primary standard candles. Regarding the SFH, the key 
contributions of the VMC survey will be: (1) It will provide 
photometry reaching as deep as the oldest main sequence turn- 
off over the bulk of the MC system, as opposed to the tiny re- 
gions sampled by HST, and the limited area covered by most of 
the dedicated ground-based observations. (2) VMC will use the 
near-infrared YJK S passbands, hopefully reducing the errors in 
the SFH-recovery due to variable extinction across the MCs. 

On the other hand, the use of near-infrared instead of op- 
tical filters will bring along some complicating factors, like an 
higher degree of contamination of the MC photometry by fore- 
ground stars and background galaxies, and the extremely high 
noise contributed by the sky, especially in the K s band. 

Indeed, VMC will be the first near-infrared wide-area sur- 
vey to provide data suitable for the classical methods of SFH- 
recoverjQ. With the new space-based near-infrared cameras 
(the HST/WFC3 IR channel, and JWST) and ground-based 
adaptive optics facilities, observations similar to VMC ones 
will likely be available for many nearby galaxies. VMC may 
become the precursor of detailed SFH-recovery in the opening 
window of near-infrared wavelengths. Demonstrating the fea- 
sibility of VMC goals, therefore, is of more general interest. 



1 See http://www.vista.ac.uk and 
http://www.star.herts.ac.uk/~mcioni/vmc/ 

for further information. 

2 The previous attempts of [Cioni et aT, I d2006allbl) based on 
IJK S data, were based on the shallow observations from DENIS 
and 2MASS, which are limited to the upper RGB and above. 
Consequently, they could access the general trends in the mean age 
and metallicity across the MCs, but not the detailed age-resolved SFH. 



Another particularity of the VMC survey is that, once 
started, its data flow will be so huge that algorithms of anal- 
ysis have better to be prepared in advance, in the form of semi- 
automated pipelines. Similar approaches have been/are being 
followed by some ambitious nearly-all-sky (SDSS, 2MASS), 
micro-lensing (MACHO, OGLE, EROS), and space astrome- 
try (e.g. Hipparcos, GAIA) surveys. 

In this paper, we describe part of the preparatory work for 
the derivation of the SFH from VMC data, which can be sum- 
marised in the following way: First we simulate the VMC im- 
ages for the LMC (Sect. |2), where we later perform the pho- 
tometry and artificial stars tests (Sect. [3) that allow us to ac- 
cess the expected levels of photometric errors, completeness, 
and crowding, and the contamination by foreground stars and 
background galaxies. We then proceed with many experiments 
of SFH-recovery (Sect. |4), evaluating the uncertainties in the 
derivation of the SFH as a function of basic quantities such as 
the stellar density over the LMC, the area included in the anal- 
ysis, and the adopted values for the distance and reddening. 
Doing this, we are able to present the expected random and 
systematic errors in the space-resolved SFH. Such information 
may be useful to plan complementary observations and surveys 
of the LMC in the next few years. Further papers will present 
the perspectives for the Small Magellanic Cloud (SMC), as 
well as better explore the effect in the recovered SFH due to 
the uncertainties associated with the MC geometry, differential 
reddening, initial mass function, fraction of binaries, etc. 

2. Simulating VMC data 

Our initial goal is to obtain realistic simulations of VMC im- 
ages, containing all of the objects that are known to be present 
towards the MCs and likely to be detectable within the survey 
depth limits. These objects are essentially stars belonging to 
the MW and the MCs, and background galaxies. Moreover, an 
essential component of the images is the high signal from the 
infrared sky. Each one of these components will be detailed be- 
low. Diffuse objects such as emission nebulae and star clusters 
will, for the moment, be ignored. 

2.1. VISTA and VMC specifications 

VMC will be performed with the VIRCAM camera mounted 
at the 4m VISTA telescope at ESO's Paranal Observatory in 
Chile. VIRCAM has 16 2048x2048 detectors which, with the 
image scale of 0.339 " per pixel on average, cover a sky area 
of 0.037 deg 2 each. The basic mode of the observations will be 
to perform 6 exposures (paw-prints) with the subsequent con- 
struction of 1 .0 x 1 .5 deg 2 tiles. In the following, we will adopt 
the area of each detector (i.e. 0.037 deg 2 ) as the basic unit of 
our simulations. 

The specifications of the VMC survey will be detailed in 
another paper (Cioni et al., in preparation). For our aims, suf- 
fice it to mention the following: Despite for the crowded fields, 
it is expected that the observations will be sky-noise dominated. 
The mean sky brightness at Cerro Paranal is of 17.2, 16.0, 13.0 
mag arcsec~ 2 in YJK S , respectively. The required seeing is of 
1 .0 " (FWHM) in the Y band, being the most crowded regions 
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Fig. 1. The area VMC will likely cover in the LMC (solid line) 
and SMC (dashed line) as a function of the surface density of 
RGB stars, /Vrgb- The fraction of the covered area in each MC, 
for four ranges of density, is also shown in the top of the figure. 
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Fig. 2. A series of Mari go et al.l (120081) isochrones in the 



UKIDSS photometric system. The figure shows the absolute 
(M K , Y-K) CMD as well as the apparent (K, Y-K) one for a 
typical distance to the LMC. Stellar masses and isochrone ages 
and metallicities are also indicated in the figure. 



observed in nights with seeing better than 0.8 ". The target- 
ted signal-to-noise ratio (SNR) is equal to 10 at magnitudes of 
21.9, 21.4, 20.3 mag respectively. The photometric zero-points 
in our simulations were fixed via the VISTA exposure time cal- 
culator, so as to be consistent with these values. Considering 
these survey limits, in our simulations we include all objects 
brighter than K s = 22.5, which at the LMC distance correspond 
to a stellar mass of ~ 0.8 M Q in the main sequence turn-off. 

VMC tiles will cover most of the Magellanic System , sum- 
ming to a total area of 1 84 deg 2 (see ICioni et all 2007 , Cioni 
et al., in preparation, for details). Fig. Q] shows a histogram 
of the total area to be observed as a function of the den- 
sity of upper RGB stars, /Vrgb, which is defined as the num- 
ber of 2MASS stars inside a box in the K s vs. J - K s CMD 
(0.60 < / - K s < 1.20 and 12.00 < K s < 14.00 for the LMC 
and 12.30 < K s < 14.30 for the SMC), for each unit area of 
0.05 deg 2 . Notice the higher mean and maximum stellar densi- 
ties of the LMC, as compared to the SMC. The stellar densities 
vary over an interval of about 2.5 dex. 



2.2. Stars in the UKIDSS photometric system 

Since VISTA is still being commissioned at the time of this 
writing, the throughputs of VISTA filters, camera, and tele- 
scope are still not available. It is however clear that the VISTA 
photometric system will be very similar to the UKIDSS one, 
with the differences being mainly in the higher performance of 
VISTA, and in the fact that VISTA will use a ^-short filter (K s ) 
similar to the 2MASS one. 



Given the present situation, we have so far used the 
UKIDSS system as a surrogate of the future VISTA one. Tests 
using the preliminary VISTA filter curves (Jim Emerson, pri- 
vate communication) indicate very small differences in the syn- 
thetic photometry, typically smaller than 0.02 mag, between 
VISTA and UKIDSS.0 

Stellar isochron es in the UK I DSS s ystem have been re- 
cently provided by IMarigo et al. ( 200 8 J. The conversion to 



the UKIDSS system takes into account not only the photo- 
spheric emission from stars, but also the reprocessing of their 
radiation by dusty s hells in mass-losing stars, as described in 
Marigo et alj ( 12008b . The filt er transmission curv es and zero- 



point definitions come from Hewett et al. (2006). The stellar 



models in use are composed of iGirardi et al. 



2000) tracks 



for low- and intermediate-mass stars, replacing the thermally 
pulsing asymptotic gian t branch (AGB) evolution with the 
Marigo & Girardil (12007b ones. In this paper, these models are 
further co mplemented with w hite and brown dwarfs as de- 
d2005l a lso Z abot et al., in prepara- 



scribed in Girardi et al 



tion), and with the Bert elli et al.l (11994 ) isochron es for masses 
higher than 7 M . Fig. [2] presents some of the IMarigo et al. 
(2008) isochrones in the Mk vs. Y-K diagram, for a wide 



range in age and metallicity. As can be readily noticed, the 



3 Throughout this paper, we will name the 2.2 fim filter as A" s when 
referring to VISTA, 2MASS and DENIS, and as K when referring 
to our simulations and to UKIDSS data. Notice however that, for all 
practical purposes, the actual difference between these filters is not a 
matter of concern. 

4 http://stev.oapd.inaf.it/dustyAGB07 
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isochrones contain the vast majority of the single objects that 
can be prominent in the near-infrared observations of the LMC, 
going from the lower Main-Sequence (MS) stars, up to the 
brightest AGB stars and red supergiants. The stellar masses 
in the MS and the appare nt magnitude for a typical LMC dis 



Grocholski et al.1 120061: iKerber et all l2007h and field stars 
dCole et al.l 120051: ICarrera et alj l2008h . together with a con- 
stant SFR(f ). Since the SFR(f) in the LMC is clearly spatial 



dependent ( Holtzma n et al. 1999; Smecker-Hane et al. 2002; 



tance, (m-M) = 18.50 (Clementi ni et alJl2003b lAlve s 2004 



ISchaeferll2008h are also indicated in this figure. We recall that 
our models contain, in addition, the very-low mass stars, brown 
dwarfs and white dwarfs, which are i mportant in the des crip- 
tion of the foreground MW population (IMarigo et al ]l2003h . 

The interstellar extinction coefficien ts adopted in this work 
do also follow from Mari go et al.l d2008l) : A 
0.283 A y , and A K = 0.114 A v , which imply E Y -j 
and Ey-K - 0.931 Eb-v- They have been derived from synthetic 



0.385 A v , Aj = 
0.351 Eb-v 



Ijaviel et al1l2005l) . the assumption of a constant SFR(f) shall 
be considered as just a way to ensure a uniform treatment for 
all stellar populations over the LMC. 

In terms of distance we are initially using the canonical 
value of (m-M )n = 18.50 dciemenfini et al.ll2003l:lAlvesll2004 : 
ISchaeferj|2008l) also adopt ed by the HST Key Pr oject to mea- 
sure the Hubble constant dFreedman et alJuOOll) . whereas for 
the reddening we are assuming a value of Eb-v = .07, typi 



photom etry applied to a G2V star extincted with lCardelli et al 



19891) extinction curve. Alth ough the approach is not the most 



cal for the extinction maps from lSchlegel et al.l (1998). For real 
VMC images, these two parameters are expected to be free pa- 
rameters, since the LMC presents disk -like geometries with a 
significant inclination (~ 30-40 deg^ dvan der Marel & Cionil 



accurate one (see |Girardi et al.|[2008|), it is appropriate to the bOOlllvan der Marel et alJ2 002: N kolaev et al. 2004) an dnon 



conditions of moderate reddening (Eb-v ~ 0.2 mag) typical of 
the Magellanic Clouds. 

The simulation of the input photometric catalogues and the 
generation of artificial images are described in the next sub- 
sections. In brief, the input catalogues for the LMC (Sect. l2~3l 
and the foreground MW stars (Sect. 12.4b come from the pre - 
dictions made with the TRILEGAL code dGirardi et al.ll2005l) . 
that simulates the photometry of resolved stellar populations 
following a given distribution of initial masses, ages, metal- 
licities, reddenings and distances. The catalogues of back- 
ground galaxies (Sect. \2.5i are randomly drawn from UKIDSS 
dLawrence et al ]l2007l) . The simulation of images is performed 
with the DAOPHOT and ARTDATA packages in IRAlf] 
(Sect. 12.6b . always respecting the photometric calibration and 
expected image quality required by the VMC survey. 

2.3. The LMC stars 

The stellar populations for the LMC are si mulated as an "addi - 
tional object" inside the TRILEGAL code dGirardi et al.l2005l) . 
where the input parameters for a field are: 

- the star formation rate as a function of stellar age, SFR(f); 

- the stellar AMR, Z(f) or [M/H](f); 

- the total stellar mass, M tot LMc; 

- the distance modulus, (ra-M)o; 

- the reddening, Eb-v = 3.1 Ay; 

- the Initial Mass Function (IMF), i^(Mj); 

- the fraction of detached unresolved binaries, /bin- 

As commented before, for convenience we are simulating 
an area of 0.037 deg 2 , equivalent to a 2048 x 2048 VIRCAM 
detector. The value of M tot ,LMC is suitably chosen such as to 
generate the total number of RGB stars observed by 2MASS, 
Nrgb, inside this same area. 

In the LMC simulations presented in this paper, we 
adopt an input AM R consistent with the one given by stel- 
lar clusters ( Olszewski et al.l [l99 lb iMackev & Gilmore 2003; 



5 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



uniformextinction dZaritskv et al. 20041 : Sub ramaniam 2005 
Imara & Blitzll2007l) . 

Finally the ass umed values for the remaining inputs are: 
the lChabrierld200ll) lognormal IMF0 , and / bin = 30 % with 
a constant mass ratio distribution for m-ilm\ > 0.7 0. There 
are no strong reasons to expect significant deviations for these 
choices, especially for the IMF since we are dealing with stars 
with masses approximately between 0.8 and 12.0 M & where the 
IMF slope seems to b e universal and similar to the Salpeter one 
dKroupall200li r2002). Concerning the fraction of binaries, our 
choice is consis tent with the values found for th e stellar clus- 
ters in the LMC dElson et al.lll998fclHu et ai1l2008l) . For the mo- 
ment, these will be considered as fixed inputs. Further papers 
will use simulations in order to quantify the systematic errors 
in the recovered SFH introduced by the uncertainties related to 
these choices. 



2.4. The Milky Way foreground 

The MW foreground stars are simu l ated u sing the TRILEGAL 
code as described in iGirardi et al.l (120051) . Towards the MCs, 
the simulated stars are located both in a disk which scale-height 
increasing with age, and in a oblate halo component. Diffuse in- 
terstellar reddening within 100 pc of the Galactic Plane is also 
considered, although it affects little the near-infrared photome- 
try. 

it has been shown that for off- 



In Girardi et al 



(12005b . 

plane line-of-sights, TRILEGAL predicts star counts accurate 
to within about 15% over a wide range of magnitudes and 
down to J a 20.5 and K ^ 18.5 . This accuracy is confir med 
by the K ^ 20.5 observations of iGullieuszik et al.l d2008l) for 



a field next t o the L eo II dwarf spheroidal galaxy. Moreover, 



Marig o et al 



d2003h shows that TRILEGAL describes very 
well the position of the three "vertical fingers" observed in 
2MASS K vs. J - K s diagrams. Similarly-comforting compar- 



-3.0 for 



6 With a slope or 2.3 for 0.8 < m/M e < 5.0 and a 

m > 5.0M o , where the Salpeter slope is a = -2.35. 

7 This is the mass ratio interval in which the secondary significantly 
affects the photometry of the system. 
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isons with UKIDSS data (including the Y band) are presented 
in Sect. |3.2| below. 

Although predicting star counts with an accuracy of about 
15% may be good enough for our initial purposes, we are work- 
ing to further improve this accuracy: In sh ort, we are applying 



the m inimisation algorithm described in IVanhollebeke et al 
(2008) - which was successfully applied to the derivation of 
Bulge parameters using data for inner MW regions - to recal- 
ibrate the TRILEGAL disk and halo models. It is likely that 
before VMC starts, foreground star counts will be predicted 
with accuracies of the order of 5%. 




2.5. The background galaxies 

In order to simulate the population of galaxies background to 
the MCs, we make use of the large catalogues of real galax- 
ies ob tained by the UKIDSS Ultra-D eep (UDS:lFoucaud et al] 
120071) and Large Area Surveys rLAS: l\Varren et alJ2007h . from 
their Data Release 3 (December 2007). The LAS includes data 
for an area of 4000 deg 2 down to K = 18.4, for YJK filters, 
whereas the UDS includes an area of 0.77 deg 2 observed down 
to K ~ 23, but only for JHK passbands. 

In our input catalogue for each simulation, we include the 
number of UKIDSS galaxies expected for our total simulated 
area. More precisely, we randomly pick up from the UDS cat- 
alogue, a fraction of galaxies given by the ratio between the 
areas covered by UDS and by our image simulation. From the 
catalogue, we extract their J and K magnitudes, and morpho- 
logical parameters (position angle, size, and axial ratio). In this 
way, our simulations respect the observed K-bwd luminosity 
function of galaxies, and their J-K colour distribution, down to 
faint magnitudes. The T-band magnitudes, instead, have been 
assigned in the following way: we take the J-K colour of each 
galaxy in the UDS, and then randomly select a galaxy from the 
LAS which has the most similar J-K (within 0.2 mag), and 
take its Y - J. This means that the Y - J vs. J - K relation from 
LAS is being extrapolated down to deeper magnitude^. 

2.6. Simulating images 

Once defined the input catalogues for stars and galaxies we 
simulated the images inside IRAF, in accordance with the 
VISTA and VMC specifications (see Sect l2.U . The basic se- 
quence of steps (and the IRAF task) performed for a given filter 
is the following: 

1 . Definition of the image size (rtextimage) and introduction 
of the sky brightness and noise (mknoise); 

2. simulation of a Gaussian stellar profile (gauss) respecting 
the expected seeing for an image of a photometric cali- 
brated (using the VISTA ETC vl.3) delta function with a 
known number of electrons; 

3. addition of the LMC and MW stars in the sky images 
(addstars) following the previous calibrated Gaussian stel- 




8 This is of course a crude approximation since deeper surveys sam- 
ple larger galaxy redshifts. However, it is justified by the lack of deep- 
enough Y data, and by the little impact that such faint galaxies have in 
our stellar photometry (see Sect. 13. 21 ). 



Fig. 3. An image simulation for the area next to the star 
cluster NGC1805 (a = 5.03 h, 5 = -66.07°), for a sin- 
gle 2048 x 2048 array detector of VIRCAM. This is a false- 
color image where blue-green-red colours were associated to 
the YJK filters, respectively. The location corresponds to a 
logA^RGB ~ 2.00 in Fig. Q] The detector area corresponds to 
0.0372 deg 2 (11.6 x 11.6 arcmin) in the sky, which is about 
1/40 of a single VIRCAM tile, and 1/5000 of the total VMC 
survey area. The two small panels at the bottom present details 
of the simulated stars, stellar clusters and galaxies for 2.9 x 2.9 
arcmin and 0.7 x 0.7 arcmin areas. At the LMC distance the 
top panel corresponds approximately to a box of 175 x 175 pc, 
whereas the bottom panels correspond to 44 x44 pc and 11x11 
pc, respectively (1 " ~ 0.25 pc). 

lar profile with random poissonian errors in the number of 
electrons; 

4. addition of galaxies (mkobject) in the previous image 
respecting all information concerning the morphological 
type, position angle, size and axial ratio. 

To assure a uniform distribution of the objects in the image, 
stars and galaxies are always added at random positions. Fig. [3] 
shows an example of image simulation for a typical field in the 
LMC disk. The false-colour plot evidences the colour and mor- 
phologic differences between stars and background galaxies - 
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with the latter being significantly redder than the former. In the 
same image, we have inserted two populous stellar clusters typ- 
ical for the LMC with different ages, masses and concentration 
of stars (following a King's profile), just to illustrate our capac- 
ity to simulate also this kind of stellar object. 

3. Performing photometry on simulated data 

3.1. Aperture and PSF photometry 

The IRAF DAOPHOT package was used to detect and to per- 
form aperture photometry in our simulated images. Candidate 
stars were detected using daofind, with a peak intensity thresh- 
old for detection set to 4 <r s k y , where cr s ky corresponds to the 
rms fluctuation in the sky counts. The aperture photometry was 
carried out running the task phot for an aperture radius of 3 pix 
(~ 1.0"). 

The photometric errors and completeness curves that come 
from this aperture photometry in our simulated LMC images 
can be seen in Fig. |U The photometric errors in this case are 
estimated using the differences between the input and output 
magnitudes; more specifically, for each small magnitude bin we 
compute the half-width of the error distribution, with respect 
to the median, that comprises 70% of the recovered stars. The 
completeness instead is simply defined as the ratio between 
total number of input stars, and those recovered by the pho- 
tometry package. Fig. |4]presents the results for different sim- 
ulations covering a wide range of density of field stars in the 
LMC, from the outer disk regions to the centre regions in the 
bar (see FigQ]). Notice that in these simulations we are follow- 
ing the requirement that the most central and crowded regions 
(log A^rgb ^ 2.50) will be observed under excellent seeing con- 
ditions only. 

It can be noticed that the VMC expected magnitudes at 
SNR=10 for isolated stars (Y = 21.9, J = 21.4, K = 20.3) 
is well recovered in the simulations for the lowest density re- 
gions, attesting the correct photometric calibration of our sim- 
ulated images. For these regions the 50% completeness level is 
reached at Y ~ 22.5, J ~ 22.2 and K ~ 21.8. 

Crowding significantly affects in the quality of the aperture 
photometry, making the stars measured in central LMC regions 
to appear significantly brighter, and with larger photometric er- 
rors, than in the outermost LMC regions. As shown in Fig. [4] 
crowding clearly starts to dominate the noise for fields with 
log A^rgb ^ 2.00 that correspond to about 20% (7%) of the 
total area covered in the LMC (SMC) (see Fig. [TJ. Therefore, 
PSF photometry is expected to be performed whenever crowd- 
ing will prevent a good aperture photometry over VMC images. 
The significant improvements that can be reach by a PSF pho- 
tometry are also illustrated by the thick black lines in Fig. [4] 
These results correspond to a PSF photometry applied to the 
LMC centre (log A^rgb = 2.90), where the PSF fitting and the 
photometry were done using the IRAF tasks psf and allstar. 

Figure shows an example of CMD obtained from the 
aperture photometry in a simulated field with an intermedi- 
ate level of density in the LMC. This figure reveals the ex- 
pected CMD features - and the wealth of information - that 
will become available thanks to the VMC survey: well evi- 




Fig. 5. Example of (K, Y-K) CMD from aperture photome- 
try in a simulated image for the VMC survey. The choices 
in the parameters represent a field of ~ 0.1 deg 2 with ~ 10 5 
stars (log(A^RGB) = 2.00) following a constant SFR(f) and a 
AMR typical for the LMC clusters (see details in the text). The 
colours represent the density of points in a logarithmic scale. 
The information about approximated stellar masses, ages and 
metallicities can be obtained from Fig. [2] 



dent are the AGB, red supergiants, RGB, red clump (RC), Sub- 
Giant Branch (SGB) as well as the MS, from the brightest 
and youngest stars down to the oldest turn-off point. In com- 
parison, the present-day near-infrared surveys of the MCs are 
complete only for the most evolved stars - excluding those in 
the most crowded regions, and those highly extincted by cir- 
cumstellar dust. DENIS and 2MASS, for instance, are lim- 
ited to K s 5= 14, revealing the red supergiants, AGB and 
upper R GB, and including just a t iny fraction of t he up- 
per MS dNikolaev & Weinbergll2000t Icioni et al.ll2007l) . IRSF 
(Kat o et all2007l) extends this range down to K s £ 16.6, which 
is deep enough to sample the RC and RGB bump, but not the 
SGB and the low-mass MS. 



3.2. Comparison with UKIDSS data 

Since the present work depends on simulations, it is impor- 
tant to check if they reproduce the basic characteristics of real 
data already obtained under similar conditions. UKIDSS rep- 
resents the most similar data to VMC to be available for the 
moment. Therefore, in the following we will compare a simu- 
lated UKIDSS field with the real one. 

For this exercise, we take the 0.21 deg 2 field towards 
Galactic coordinates I = -220°, b = 40°, which due to its 
similar distance from the Galactic Plane as the MCs, offers a 
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magnitude,, magnitude. 

o in o m 

Fig. 4. Photometric errors (left panels) and completeness curves (right panels) in the artificial YJK images for the LMC for 
different levels of crowding. The thin curves present the results for the aperture photometry covering the entire expected range 
of density of field stars (logjV RGB = 1-50, 1.75, 2.00, 2.25, 2.50, 2.75, 2.90, see also Fig. [D- The three highest density levels 
were simulated with the smallest values for seeing required for the LMC centre. The thick black line illustrates the results of 
performing PSF photometry for the highest density level (log A^rgb = 2.90). The expected error in magnitude for a SNR=10 is 
shown by the dashed line. 



good opportunity to test the expected levels of Milky Way fore- 
ground, and the galaxy background. 

We have taken the original image from the UKIDSS 
archive, an d performed aperture photo metry with both 



DAOP HOT dStetsonl fl987l) ad SExtractor feerfin & Arnouts 
1996h . A image for the same area has been simulated using 
UKIDSS specifications (pixel scale, SNR, etc.) and submitted 
to the same kind of catalogue extraction. Fig. [6] shows the re- 
sults, comparing the K vs. Y-K diagrams for the UKIDSS ob- 
served (left panel) and simulated (right panel) fields, for both 
stars (blue points) and galaxies (red points). Stars and galaxies 
were separated using the SExtractor Stellarity parameter st. We 
adopted st > 0.85 for stars, and st < 0.85 for galaxies. 

Both DAOPHOT ad SExtractor aperture photometries 
turned out to be remarkably consistent with the ones provided 
by the Cambridge Astronomical Survey Unit (CASU) data re- 
duction pipeline. This is very comforting since the CASU will 
adopt the same data reduction pipeline to the future VISTA 
data. The histograms at the right and top of the CMD pan- 
els show the object number count distribution in both colour 
and magnitude. As can be appreciated, our simulated objects 
distribute very similarly in colour and magnitude as the ob- 
served ones. The discrepancies are limited to a few aspects of 
the simulations, for instance the peak in the colour distribution 
at Y-K ~ 1.3 is clearly narrower in the models than in the 
simulations. This peak is caused by thin disk dwarfs less mas- 
sive than 0.4 M Q dMarigo et al.ll2003 ). and its narrowness in the 
models could be indicating that TRILEGAL underestimates the 
colour spread of these very-low mass stars. 



The most important point of the model-data comparison of 
Fig. [6] however, is that the simulations reasonably reproduce 
the numbers (with errors limited to ~ 20%), magnitudes and 
colours of the observed objects. This gives us confidence that 
our MC simulations contain the correct contribution from fore- 
ground Milky Way stars and background galaxies. 



4. Recovering the SFH 

4.1. Basics 

The basic assumption behind any method to recover the SFH 
from a composite stellar population (CSP) is that it can be 
considered as a simply sum of its constituent parts, which 
are ultimately simple stellar populations (SSPs) or combina- 
tions of them. Therefore determining the SFH of any CSP - 
like the field stars in a galaxy - means to recover the rela- 
tive weight of each SSP The modern stellar population analy- 
sis in the late 8 0's dTosi et al.ll 1 989b iFerraro et all 1989b and in 
the early 90's dTosi et alJll99lb iBertelli et al.lll992h - marked 
by the advent of the first CCD detectors - was done assum- 
ing parameterized SFH, which revealed the main trends in the 
SFH but was still limited by a small number of possible solu- 
tions. The techniques to recover the SFH from a resolved stel- 
lar popul ation started t o become more sophisticated with the 
works of Galla rt et al. I dl996bl aj), b ut t hey were significantly 
improved by lAparicio et al.l d 19971) and Dolphinl (11997b . who 
developed for the first time statistical methods to recover non- 
parameterized SFH from the CMD of a CSP. In practice these 
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Fig. 6. A comparison between the aperture photometry from UKIDSS image data (left) and the corresponding simulation (right), 
for a 0.21 deg 2 area towards I = -220°, b = 40°. The photometry was performed using both DAOPHOT and SExtractor. The 
main panels show the CMD obtained combining DAOPHOT photometry with SExtractor star/galaxy classification (blue/red 
dots, respectively). The histograms show the total colour and magnitude distributions of stars and galaxies (blue and red lines, 
respectively). 



two works were the first to deal with a finite number of free 
independent components, obtained by adding the properties of 
SSPs inside sma ll, but finite, age and metallicity bins. These 
"partial models" (lAparicio et alJ 19971) are so computed for age 



and metallicity bins that should be small enough so that the SSP 
properties change just little inside them, and large enough so 
that the limited number of bins ensures reasonable CPU times 
for the SFH-recovery. Furthermore, being the partial models 
computed for the same and constant star formation rate in- 
side each age bin, implied that they needed to be generated 
only once, sav ing a large amount of computational resources 
dDolphinl2002h . 

Considering that a CMD is a distribution of points in a 
plane divided into AW boxes, these ideas can be expressed 
ilDolphinl2002l) by 



= Z r i Ci - 



(1) 



where m, is the number of stars in the full model CMD for an 
CSP in the / th box, r> is the SFR for the / h partial model, and 
Cu is the number of stars in the CMD for the / h partial model 
in the i th CMD box. 

Notice that the above equation is in fact written in terms of 
Hess diagrams since we are dealing with the number of stars 
in CMDs. So it means that the "observed" Hess diagram for a 
CSP can be described as a sum of independent synthetic Hess 
diagrams of partial models, being the coefficients r s the SFRs to 
be determined. Fig. [71 to be commented later, illustrate the gen- 
eration of such synthetic Hess diagrams for the partial models 
of the LMC. 



The classical approach to determine the set of rjs is 
to compute the differences in the number of stars in 
each CMD box between data and model, searching for 
a minimisation of a chi-squared-like statistics. This kind 
of approach was applied f or the first time to r ecover the 
SFH of a real galaxy by lAparicio et al.l dl997l) . and has 
been successfully used in the analysis of th e field stars 
in th e dwarf galaxie s in the Loca l Grou p dGallart et al 



1 1999 1 : iDotohiri |2002[ fPolphin et al] 120031: ISkillman et a" 
120031: ICole elaLlboOTt l?uk & Leell2007l). including the MCs 



(Olsed [l999UHoltzman et al. 19991 lHarris & Zaritsk\] 2001 



Smecker-Hane et alj200~ ; Harris & Zaritskvl2004l:|javiel et al 

alJ| 



120051 IChiosi & Vallenaril 120071: iNoel et alJ 120071) . Although 
these works have in common the same basic idea of how to re- 
cover the SFH, there are clear variations in terms o f the adopted 



statis tics and strategy to divide the CMD in boxes dGallart et al 
l2005h . 

An interesting alternative to recover the SFH from the 
analysis of CMDs is offered by the maximum likelihood 
technique using a Bayesian approach dTolstov &Sahalll996i: 
Hernandez etaLlll999L Eoool IVergely et al.ll2002b . In this ap- 



proach the basic idea is to establish for each observed star 
a probability to belong to a SSP (based on the expected 
number of stars from this SSP in the position of the ob- 
served star in the CMD). Doing it for all observed stars, 
one can recover the SFRs that maximises the likelihood be- 
tween data and model. It is interesting to note that in the re- 
cent years there is an increasing number of papers applying 
this kind of technique for a wide range of problems, which 
include the determination of physical parameters of stellar 
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Fig. 7. Simulated {K, Y-K) CMD illustrating the building of partial models for the analysis of LMC stellar populations and its 
foreground MW stars. Panel (a) shows the theoretical stars generated from TRILEGAL, corresponds to the following ranges in 
log(f/yr): 8.00-8.40 (red), 8.40-8.80 (green), 8.80-9.20 (blue), 9.20-9.60 (cyan), 9.60-10.00 (magenta), 10.00-10.15 (yellow), 
plus the foreground MW (black). Panel (b) shows the same after considering the effects of photometric errors and completeness. 
Panel (c) is the Hess diagram for the sum of all partial models. P anels (d) and (e) show the input SFR(f) and AMR, respectively, 
the latter in comparison with LMC clusters (squares - data from Mackey & Gilmorel l2003 : iKerber et aL 2007 : Grocholski et al.l 
2006, 2007). To avoid an extremely large size for this figure only 5% of all stars typically used (~ 10 7 ) to build the partial models 
are shown in the top panels. 



clusters dfergensen & Lindegrenl2005 : Nay lor & Jeffriesl2 006: 
Hernandez & Valls-Gabaud 2008) as well as of individual stars 
dNordstrom et al.ll2004t Ida Silva et ai1l2006l) . 

It is beyond the scope of the present work to discuss in 
depth the particularity of each aforementioned approach, but 
there are no strong reasons to believe that o ne can intrinsi- 
cally recover a mor e reliable SFH than the other ( Dolphin 2002; 
Gallart et al.ll20 05). So for a question of simplicity and coher- 
ence with the majority of the works devoted to the MCs, we 
adopted the classical minimisation of a chi-squared-like statis- 
tics technique to determine the expected errors in the SFH 



for the VMC data, using the fra mework of the StarFISH code 
dHarris & Zaritskyll200ll |2004 . the ^ 2 -like statistics defined 
bv lDolphinl 020021) assuming that stars into CMD boxes follow 
a Poisson-distributed data, and a uniform grid of boxes in the 
CMD. 



4.2. StarFISH and TRILEGAL working together 

The S tarFISH cod|] has been devel oped by lHarris & Za ritsk' 



ine M arriaH coaqj has been developed by Harris 8i z,aritskv 
( 2001 ) and successfully applied by Harris & Zaritskvl "( 20041) 



Available at http://www.noao.edu/staff/jharris/SFH/ 
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and Harris (2007b) to recover SFHs for the MCs inside 
the context of the Magellanic Clouds Photometric Survey 
dZaritskvet al.ll 19971 MCPS0 This code, originally designed 
to analyse CMDs built w ith UBVI data from the M CPS, and 
using Padova isochrones dGirardi et al]|2000l 2002b . offers the 
possibility of different choices for generating synthetic Hess di- 
agrams (set of partial models, CMD binning and masks, com- 
bination of more than one CMD, etc.) and^ 2 -like statistics, be- 
ing also sufficiently generic to be implemented for new stellar 
evolutionary models, photometric systems, etc. 

As illustrated in Fig. [7] the TR1LEGAL code can also sim- 
ulate the synthetic Hess diagram for a set of partial models, 
with the advantage of easily generating them in the UKIDSS 
photometric system, allowing also a greater control of all input 
parameters involved (see Sect. 12.3b . Therefore we decide to di- 
rectly provide these Hess diagrams to StarFISH, using it as a 
platform to determine the SFRs for our VMC simulated data 
by means of a ^ 2 -like statistics minimisation. The search of 
the best solution is done internally in StarFISH by the amoeba 
algorithm that uses a downhill strategy to find the minimum 
^ 2 -like statistics value. 

An extra possibility offered by TRILEGAL is the construc- 
tion of an additional partial model for the MW foregroun 

dB 

Indeed, this is done by simulating the MW population towards 
the galactic coordinates under examination, for the same to- 
tal observed area but averaging over many simulations so as 
to reduce the Poisson noise. This partial model is provided to 
StarFISH and used in the ^ 2 -like statistics minimization to- 
gether with those used to describe the MC population. With this 
procedure, the presence of the MW foreground is taken into ac- 
count in the SFH determination, without appealing to the (often 
risky) procedures of statistical decontamination based on the 
observation of external control fields. To our knowledgement, 
this is the first time that such a procedure is adopted in SFH- 
recovery work. Notice that, once the MW foreground model 
is well calibrated, its corresponding rj could be set to a fixed 
value, instead of being included into the^ 2 -like statistics min- 
imization. 

Figure|7]illustrates the generation of a complete set of par- 
tial models, covering ages from log(f/yr) = 8.00 to 10.15 (? 
from 0.10 to 14.13 Gyr) divided into 1 1 elements with a width 
of A log t = 0.20 each, and following an AMR consistent with 
the LMC clusters (see the panel d, and Sect. 12.31 ). In this figure 
we have grouped the partial models in just 6 age ranges (plus 
the MW foreground one) just for a question of clarity. What is 
remarkable in the figure is the high degree of superposition of 
the different partial models over the RGB region of the CMD 
- except of course for the partial model corresponding to the 
MW foreground. The MS region of the CMD, instead, allows a 
good visual separation of the different populations over the en- 
tire age range, even after considering the effects of photometric 
errors and incompleteness. 
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Fig. 9. Errors in the recovered SFR(f) in terms of the mean 
SFR(f) (top panel) and input SFR(f) (bottom panel). The in- 
put simulations correspond to a typical LMC disk region 
(log /Vrgb = 2.00) inside a single VIRCAM detector (~ 
0.037 deg 2 ). The central solid line corresponds to the median 
solution found over 100 realizations of the same simulation, 
whereas the error bars correspond to a confidence level of 70%. 



4.3. Results: Input vs. output SFR{t) 

An example of SFH-recovery is presented in the Hess diagrams 
of Fig. [8] The input simulation (left panel) was generated for a 
constant SFR(f), for an area equivalent to 1 VIRCAM detector 
(0.037 deg 2 ) inside a region with a stellar density typical for 
the LMC disk (log /Vrgb = 2.00, that produces a total number 
of stars of /V stars ~ 5 x 10 4 ). For this simulation, StarFISH fits 
the solution represented in the middle panel. Not surprisingly, 
the data-model ^- 2 -like statistics residuals (right panel) are re- 
markably evenly distributed across the Hess diagram. 

Figure [9] presents the median of recovered SFR(f), and its 
error, obtained after performing 100 realizations of the same 
simulation. As expected, the median SFR(f) over this 100 real- 
izations reproduces remarkably well the input one, with no in- 
dication of systematic errors in the process of SFH-recover){3 
. The error bars correspond to a confidence level of 70%, 
which means that 70% of all individual realizations are con- 
fined within these error bars. Error bars are almost symmetri- 
cal with respect to the expected SFR(t). Furthermore, errors are 
typically below 0.4 in units of mean SFR(f) (top panel), which 
means uncertainties below 40% (bottom panel) for almost all 
ages. The only exception is the youngest age bin which presents 
errors in the SFR that are about two times larger than those for 
the intermediate-age stellar populations. 



10 http ://ngala. as . arizona. edu/dennis/mcsurvey . html 

11 See http://stev.oapd.inaf.it/trilegal 



12 As a consequence, the total integral of the SFR is also well recov- 
ered. 
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Fig. 8. Example of simulated 
(K, Y-K) Hess diagrams for the 
VMC (left panel) and for the best 
solution found by StarFISH (mid- 
dle panel). The ;^ 2 -like statistics 
map is also shown in the right 
panel. These Hess diagrams are 
limited to 12.00 < K < 20.50 and 
-0.50 < Y-K < 2.20 and were 
built with bin sizes of 0.10 mag 
both in colour and in magnitude. 



There are a many factors that can affect the accuracy of a 
recovered SFR(f). Among them, the most important are: 

1 . the quality of the data in terms of stellar statistics and pho- 
tometry, that depends in principle (for the same photomet- 
ric conditions of seeing, exposure times, calibration, etc.) 
on the density of the field and its covered area in the sky; 

2. the uncertainties in the models themselves, that come from 
the possible errors in the adopted input parameters (dis- 
tance, reddening, IMF, /bin, AMR, etc.), in the stellar evo- 
lutionary models, and in the imperfect reproduction of pho- 
tometric errors and completeness; 

3. the incorrect representation of the contamination from 
other sources, like foreground MW stars, stars from LMC 
star clusters, and background galaxies; 

4. the non-uniform properties of the analysed field, like differ- 
ential reddening or depth in distance. 

Notice that the first factor affects the generation of observed 
Hess diagrams, while the synthetic Hess diagrams may become 
unrealistic due to the other factors. Furthermore, they produce 
different types of errors: whereas the first preferentially rules 
the random errors, the second is the main source of the system- 
atic errors. A discussion on errors is addressed below. 

4.3.1 . Random errors for a known AMR 

In order to estimate the expected random errors in the SFR(f) 
for the LMC we performed controlled experiments similar to 
the one shown in Fig. [9] but covering a wide range of condi- 
tions in regard to the stellar statistics and crowding. Fig. [10] 
presents the results for four different levels of density, from the 
outer LMC disk (top-left panel) to the LMC centre (bottom- 
right panel), for a number of stars (or area) varying by a factor 
of 8. As can be seen, these two factors dramatically change the 
level of accuracy than can be achieved in the recovered SFR(f): 
an increase in the number of stars reduces the errors whereas 
an increase in crowding for a fixed number of stars acts in the 
opposite way. 




9.5 10 8 8.5 

log(r/yr) 

Fig. 10. Errors in the recovered SFR(f) for four different stel- 
lar densities, from the outer LMC disk (top-left panel) to the 
LMC centre (bottom-right panel). The thick black solid line 
corresponds to the median solution found using 4 x 10 5 stars, 
whereas the error bars correspond to a confidence level of 70%. 
The thin red lines outline the same confidence level for a de- 
creasing number of stars: 2 x 10 5 , 10 5 , 5 x 10 4 stars. 

The errors in the recovered SFR(f) as a function of the cov- 
ered area, for all simulated density levels, are shown in Fig.fTTI 
for partial models of four different ages, from young (top-left 
panel) to old ones (bottom-right panel). Here, a very interest- 
ing result can be seen: for stars older than log(f/yr) ~ 8.60 
(f ~ 0.4 Gyr) the curves for different levels of density are al- 
most superposed, revealing that for a fixed area the accuracy 
in the recovered SFR(f) is roughly independent of the level of 
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Fig. 11. Errors in the recovered SFR(f) as a function of covered 
area and density of stars (different lines) for 4 partial models of 
different ages (different panels). 

crowding. It can be understood as a counterbalanced effect be- 
tween the loss of stars due to a decrease in completeness, and 
the gain of stars due to an increase in density. Therefore, it 
seems that the SFR(f) for this wide range in age can be deter- 
mined with random errors below 20% if an area of 0.10 deg 2 
is used. Increasing the area by a factor of four means that the 
level of uncertainty drops to below 10%. 

On the other hand, for partial models younger than 
log(f/yr) ~ 8.60 (f ~ 0.4 Gyr) a fixed area does not imply a 
constant level of accuracy in the SFR(f), being the errors sig- 
nificantly greater in the less dense regions. Since stars in this 
small age range are mainly identified in the upper main se- 
quence at 18 ^ K % 20, and in the core-helium burning phases 
at 14 S5 K S= 18 (see Fig. [7]), this effect can be understood by the 
fact that for these brighter stars the increase in stellar density is 
not followed by a significant decrease in completeness. Indeed, 
Fig.|4]indicates the completeness is in general higher than 60% 
and 90%, for 18 < K < 20 and 14 < K < 18 stars, respectively, 
for the entire range of stellar densities of the LMC. In this sit- 
uation, what determines the accuracy for a fixed area is sim- 
ply the number of observed stars, which is proportional to the 
density. In particular, our simulations reveal the lack of stellar 
statistics in the outermost and less dense LMC regions, which 
require great areas to reach a statistically significant number of 
young stars. For instance, a level of uncertainty 20 % in the re- 
covered young SFR(f) is obtained for an area ~ 10 times larger 
(~ 1 deg 2 ) in the periphery of the LMC than in more central 
regions. 

The errors in the SFR(f) were also compared with those 
expected by the Poisson statistics in the number of observed 
stars for each partial model, as presented by Fig. [12] This figure 
reveals that the errors in the SFR(f) can be roughly understood 



Fig. 12. Distribution of the fraction of observed stars (top 
panel) and its Poisson relative errors (bottom panel, black lines) 
as a function age. These errors correspond to 4 x 10 5 simu- 
lated stars. The relative errors in the recovered SFR(f) are also 
shown (bottom panel, red lines) for a vertical scale 10 times 
larger (compare the labels of the two vertical axis in the bottom 
panel). 



as a propagation of the Poisson errors by a factor -10. This 
explains why errors are smaller for the older ages, despite the 
fact that the different partial models are - due to the adoption 
of a logarithmic age scale - roughly uniformly separated in the 
CMD. 

4.3.2. Random errors for an unknown AMR 

In our previous discussion we have naively assumed that the 
AMR of the LMC is well known, by using a set of partial mod- 
els which strictly follows the AMR used in the simulations. The 
real situation is much more complicated. Not only the AMR is 
not well established, but also it may present significant spreads 
(for a single age), and vary from place to place over the LMC 
disk. In order to face this situation, it is advisable to allow a 
more flexible approach for the SFH-recovery, in which for ev- 
ery age bin we have different partial models covering a signifi- 
cant range in metallicity. 

We adopt the scheme illustrated in the top panel of Fig. [T3j 
that is: for every age bin, we build partial models of 5 different 
metallicities, centered at the [M/H] value given by the reference 
AMR, and separated by steps of 0.2 dex. This gives a total of 
56 partial models, and drastically increases the CPU time (by a 
factor of ~ 100) needed by StarFISH to converge towards the 
^ 2 -like statistics minimum. Once the minimum is found, for 
each age bin we compute the r,(f)-added SFR(f) and the rj(t)- 
weighted average [M/H](f) over the five partial models with 
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Fig. 14. Errors in the recovered SFR(f) in terms of the mean SFR(f) (top panels) and input SFR(f) (bottom panels), for a typical 
LMC disk region inside the area of 8 VIRCAM detector (0.25 deg 2 ). The left panel correspond to a SFH-recovery that uses 
partial models distributed over a single AMR, whereas the right one uses 5 partial models for each age bin. The central solid line 
corresponds to the median solution found over 100 realizations of the same simulation, whereas the error bars correspond to a 
confidence level of 70%. 



different levels of metallicity. After doing the same for 100 re- 
alizations of the input simulation, we derive the median and the 
confidence level of 70% (assumed as the error bar) for the out- 
put SFR(f) and [M/H](f). Figs. Qj] and [14] illustrate the results 
for the case of a constant input SFR(f) and log Nrgb = 2.00. 



In the top panel of Fig. Q~3] the dots with error bars show 
the output AMR, which falls remarkably close to the input one 
(continuous line). The error bars are smaller than 0.1 dex at all 
ages. The bottom panel plots the relative errors in the derived 
[M/H](f), showing that they are slightly larger for populations 
of age log(f/yr) < 9.2. Anyway, the main result here is that the 
errors in [M/H](f) are always smaller than the 0.2 dex separa- 
tion between the different partial models. 



Figure [14] instead compares, for the same simulation, the 
errors in the SFR(f) that result either considering (right panel) 
or not (left panel) the partial models with metallicity different 
from the reference AMR one. In other words, the left panel 
shows the SFR(f) that would be recovered if the AMR were 
exactly known in advance, whereas the right panel shows the 
SFR(f) for the cases in which the AMR is unknown - or, al- 
ternatively, affected by significant observational errors. It can 
be easily noticed that the SFR(f) is correctly recovered in both 
cases, although errors in the second case (right panel) are about 
2 or 3 times larger than in the first case. Needless to say, the 
second situation is the more realistic one, and will likely be the 
one applied in the analysis of VMC data. 



4.3.3. Systematic errors related to distance and 
reddening 

Accessing all the systematic errors in the derived SFR(f) is be- 
yond the scope of this paper. However, the errors associated 
with the variations in the distance and reddening are of par- 
ticular interest here, since both quantities are expected to vary 
sensibly across the LMC, hence having the potential to affect 
the patterns in the spacially-resolved SFH. Fortunately, these 
errors are very easily accessed with our method, since we know 
exactly the (m— M)o and Ay values of the simulations, as well 
as those assumed during the SFH-recovery. 

The propagation of the uncertainties in the assumed dis- 
tance modulus and reddening in the recovered SFR(f) was ex- 
plored in a series of control experiments, in which a simu- 
lation performed with the canonical (m— M)o = 18.50 and 
Eb-v = 0.07 was submitted to a series of SFH-recovery anal- 
yses covering a range of (m-M)o and Eb-v- More specifically, 
(m-M)o was varied from 18.40 to 18.60, and Eb-v from 0.04 
to 0.1. The results of these tests can be seen in Fig. [15] which 
presents the effect of wrong choices in these parameters not 
only in the recovered SFR(f) (bottom panel) but also in the min- 
imum value for the ;f Mike statistics (top panels). 

As expected, the absolute minimum value for the ^ 2 -like 
statistics was found for the synthetic Hess diagrams with the 
right distance modulus and reddening. Based on the ^ 2 -like 
statistics dispersion for 100 simulations we estimated that the 
errors in these parameters are about A(ra-M)o = ±0.02 mag 
and AEb-v = ±0.01 mag. These errors also imply systematic 
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Fig. 13. The top panel shows the distribution of metallicities 
of the partial models adopted in this work: The central solid 
line is the AMR adopted as a reference, and is used in all of 
our SFH-recovery experiments. At every age (or age bin), 4 
additional partial models (along the dotted lines) can be defined 
and inserted in the SFH-recovery, then allowing us to access 
the AMR and its uncertainty (see text for details). The bottom 
panel shows the difference between the input and output AMRs 
for the case in which 5 partial models were adopted for each 
age bin. 



errors in the recovered SFR(f) of up to ~ 30% (Fig. Q3] bottom 
panel). 

The above experiments clearly suggest us that the mean 
distance and reddening should be considered as free parame- 
ters in the analysis, and varied by a few 0.01 mag so that we 
can identify the best-fitting values of (m-M)o and Eb-v to- 
gether with the best -fitting SFH. Th i s kind of procedure has 
been adopted by e.g. lHoltzman et al.l (11997b and lOls en ( 1999) 
in their study of small regions over the LMC, and it is also im- 
plemented in the MATCH SFH-recovery package by Dolphin 
(2002). Occasionally, one could also consider small spreads in 
both (m-M)o and Eb-v and test whether they further improve 
the ^ 2 -like statistics minimization. Once applied to the entire 
VMC area, the final result of this procedure will be indepen- 
dent maps of the geometry and reddening across the MC sys- 
tem, that can complement those obtained with other methods. 



5. Concluding remarks 

In this paper we have performed detailed simulations of the 
LMC images expected from the VMC survey, and analysed 
them in terms of the expected accuracy in determining the 
space-resolved SFH. Our main conclusions so far are the fol- 
lowing: 
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Fig. 15. The top panels show the minimum ^ 2 -like statistics 
value as a function of distance modulus (blue) and reddening 
(green) adopted to build the synthetic Hess diagrams. The bot- 
tom panel illustrate the systematic variations that the errors in 
these two parameters cause in the recovered SFR(f) (continu- 
ous lines for an overestimation of (m - M)q and Eb-v, dashed 
lines for an underestimation). 



1 . For a typical 0.10 deg 2 LMC field of median stellar density, 
the random errors in the recovered SFR(f) will be typically 
smaller than 20% for 0.2 dex-wide age bins. 

2. For all ages larger than 0.4 Gyr, at increasing stellar den- 
sities the better statistics largely compensates the effects of 
increased photometric errors and decreased completeness, 
so that good-quality SFR(f) can be determined even for the 
most crowded regions in the LMC bar. The SFR(f) errors 
decrease roughly in proportion with the square of the to- 
tal number of stars. The exception to this rule regards the 
youngest stars, which because of their brightness are less 
affected by incompleteness. In this latter case, however, 
the stellar statistics is intrinsically small and large areas 
are necessary to reach the same SFR(f) accuracy as for the 
intermediate-age and old LMC stars. 

3. Although the AMR [M/H](f) can be recovered with accu- 
racies better than 0.2 dex, the uncertainties in the AMR 
can significantly affect the quality of the derived SFR(f), 
increasing their errors by a factor of about 2.5. 

4. The minimisation algorithms allow to identify the best- 
fitting reddening and distance with accuracies of the order 
of 0.02 mag in distance modulus, and 0.01 mag in Eb-v- 

All of the above trends were derived from analysis of small 
LMC areas, that we have considered to be homogeneous in 
all of their properties (AMR, distance and reddening). The er- 
rors were derived by varying each one of these parameters 
separately. The real situation will be, of course, much more 
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complicated, with significant spatial variations of all of these 
quantities across the LMC. This consideration may lead us to 
suppose that errors here derived are underestimated. However, 
the above-mentioned parameters can be further constrained by 
simply taking into consideration additional data - for instance 
the available reddening maps, the limits on the relative dis- 
tances provided by other independent distance indicators, the 
metallicity distributions of field stars, etc. - in our analysis. 
Moreover, our work indicates clearly how the random errors 
are reduced when we increase the area to be analysed. It is nat- 
ural that, once the systematic errors are fully assessed, we will 
increase the area selected for the analysis, so that random errors 
become at least smaller than the systematic ones. 

It is also worth mentioning that our present results were 
obtained using the Y—K colour only. VMC will also provide 
CMDs involving the J passband, and their use in the SFH anal- 
ysis can only reduce the final errors. 

Another factor to be considered, in the final analysis, is 
that for old ages the SFR(f) is expected to vary very smoothly 
acro ss the LMC, as indicat e d for i nstance bv lCioni et alj yOOO) 
and Nikolaev & Weinberg (2000). This large-scale correlation 
in the old SFR(f) may be used as an additional constraint dur- 
ing the SFH-recovery, and may help to reduce the errors in the 
SFR(r) at all ages. 

A forthcoming paper will discuss in more detail the ex- 
pected accuracy over the complete VMC area, including the 
SMC, and how this accuracy depends on other variables and 
constraints which were not discussed in this work (binaries, 
depth in distance, reddening variations, IMF, etc.). Anyway, 
the present work already illustrates the excellent accuracy in 
the measurements of the space-resolved SFH, that will be made 
possible by VMC data. Moreover, it demonstrates that detailed 
SFH-recovery using deep near-infrared photometry is also fea- 
sible, as much it has always been for the case of visual obser- 
vations. 
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